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SIMULATION APPARATUS AND SIMULATION METHOD 

CROSS REFERENCE TO RELATED APPLICATIONS 
This application is based upon and claims the benefit 
5 of priority from the prior Japanese Patent Applications No. 
2001-362629, filed on November 28, 2001, the entire contents 
of which are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
10 Field of the Invention 

The present invention relates to a simulation apparatus 
and a simulation method for estimating properties of a 
semiconductor device by taking into consideration a stress 
p having occurred during the fabrication of the semiconductor 

m 15 device. 

m 

*{ Related Background Art 

With the progress of micro fabrication technology, it 

p has become impossible to develop semiconductor devices 

y- 

pj 20 without taking into consideration the effects given by a 

stress occurring during fabrication of the semiconductor 
devices. For example, there is a problem that the mobility 
of electrons becomes low due to a compression stress occurring 
in a substrate, thereby decreasing the current value of an 
25 N-MOSFET as compared with the case where the compression 
stress is lower. 

On the other hand, a P-MOSFET has an advantage that the 
current value increases in accordance with the increase in 
the mobility of holes, thereby improving device properties. 
30 Besides the above-described matters, because the 

stress occurring in the semiconductor device has influence 
on the impurity diffusion and the shape of the semiconductor 
device, a semiconductor process device simulation capable 
of easily taking a countermeasure for the stress is important. 
35 In the conventional semiconductor process device 

simulation, the value of stress occurring during the 
fabrication of the semiconductor device is calculated, a 
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stress distribution in the semiconductor device is displayed, 
and based on the displayed result, effects of the stress on 
the device properties has been estimated experientially . 

For example, Fig. 9A shows a conventional example of 
5 displaying the stress distribution with contours in an area 
right under the gate shown in Fig. 9B. Fig. 9A displays the 
contours of the stress value based on the stress value 
calculated by using a stress simulator. 

Because the stress distribution shown in Fig. 9A is only 
10 information for estimating the effect of the stress, it was 
difficult to quantitatively estimate the effect of the stress 
on the device properties. 

SUMMARY OF THE INVENT IO N 
15 A simulation apparatus configured to estimate 

properties of a semiconductor device, comprising: 

a first calculating part configured to calculate a 
first value corresponding to a prescribed physical property 
value by taking a prescribed physical quantity into 
20 consideration, with regard to at least a partial region of 
said semiconductor device; 

a second calculating part configured to calculate a 
second value corresponding to said physical property value 
without taking said physical quantity into consideration, 
25 with regard to at least a partial region of said semiconductor 
device, and 

a visualizing part configured to display, in a 
prescribed form, a correlation between said first and second 
values . 

30 

B RIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram showing schematic 
configuration of a first embodiment of a simulation apparatus 
according to the present invention. 
35 Fig. 2 is a flow chart showing processes of the 

simulation apparatus of Fig. 1. 

Fig. 3A is a diagram showing the stress distribution 
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in a silicon substrate right under the gate. 

Fig. 3B is a diagram showing an example of the changing 
amount jLL{<7)//U of the mobility of the electrons. 

Fig. 4 is a flow chart showing processes of the second 
5 embodiment of a simulation apparatus according to the present 
invention. 

Fig. 5A shows the stress distribution in the silicon 
substrate right under the gate, similarly to Fig. 3A, and 
Fig. 5B is a diagram showing the process result of the step 
10 S16 of Fig. 4. 

Fig. 6 is a flow chart showing processes of the third 
embodiment of a simulation apparatus according to the third 
embodiment . 

Fig. 7A is the stress distribution occurring in the 
15 silicon substrate of the bulk, Fig. 7B is the stress 
distribution occurring in the SOI substrate, and Fig. 7C is 
a diagram showing the changing amount of the mobility of the 
electrons calculated from the stress distribution of Fig. 
7B. 

20 Fig. 8 is a diagram showing overall configuration of 

a simulation apparatus according to the present invention 

in the case of using a computer. 

Fig. 9A is a diagram showing the conventional example 

for displaying with the contour line the stress distribution 
25 right under the gate, and Fig. 9B is a top face view of the 

semiconductor device. 



DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Hereinafter, a simulation apparatus and a simulation 
30 method according to the present invention will be 
specifically described with reference to drawings. 
(First Embodiment) 

Fig. 1 is a block diagram showing schematic 
configuration of a first embodiment of a simulation apparatus 
35 according to the present invention. The simulation 
apparatus of Fig. 1 has a device simulation executing part 
1 for executing a simulation for a property estimation of 
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a semiconductor device based on a prescribed fabrication 
condition, a stress simulation executing part 2 for executing 
a stress simulation occurring in the semiconductor device, 
an input part 3 for inputting a fabrication condition and 
5 so on for the simulation, an output part 5 having a display 
apparatus 4 for outputting the simulated result, a first 
calculating part 6 for calculating a physical property value 
such as a mobility by taking into consideration the stress 
distribution in the semiconductor device, a second 

10 calculating part for calculating the physical property value 
without taking the stress distribution in the semiconductor 
device, a physical property value changing amount calculating 
part 8 for calculating the changing amount of the physical 
property value obtained by taking the stress into 

15 consideration with regard to the physical property value 
obtained without taking the stress into consideration, and 
a visualizing display control part 9 for visually displaying 
the result calculated by the physical property changing 
amount calculating part 8 . 

20 The device simulation executing part 1 executes the 

simulation with regard to the fabrication processes such as 
oxidation, diffusion, deposition, etching and ion 
implantation of the semiconductor device. 

The physical property values calculated by the first 

25 and second calculating parts 6 and 7 is, for example, a carrier 
mobility, carrier trap, a fixed electrical charge, a tunnel 
probability, a lifetime of carrier, a generation-extinction 
speed of pair of carriers, a diffusion coefficient, a viscous 
coefficient, and physical property values relating thereto. 

30 Fig. 2 is a flow chart showing processes of the 

simulation apparatus of Fig. 1. First of all, at least a 
partial region of the semiconductor device is divided on mesh 
form (step SI). Subseguently , the stress simulation 
calculation executing part 2 calculates the stress on each 

35 of grid points divided on mesh form, and calculates the stress 
distribution of the semiconductor device based on the 
calculated result of the stress (step S2 ) . 



Fig. 3A is a diagram showing the stress distribution 
in a silicon substrate right under the gate. In Fig. 3A, the 
deeper the color is, the stronger the compression stress is, 
and the lighter the color is, the stronger the tensile stress 
is. As shown in Fig. 3A, the nearer to the surface the region 
is, the stronger the compression stress is. The stress of 
Fig. 3A expresses a hydrostatic pressure obtained as an 
average between the stress in a crosswise direction and the 
stress in a lengthwise direction. 

Next to the step S2 , the physical property value for 
looking into effects due to the stress is selected (step S3) . 
Here, as described above, the carrier mobility and so on are 
selected. 

Subsequently, the first calculating part 6 calculates 
the physical property value selected by the step S3 by taking 
the stress distribution into consideration. Furthermore, 
the second calculating part 7 calculates the physical 
property value selected by the step S3 without taking the 
stress distribution into consideration step S4). 

Subsequently, the physical property value changing 
amount calculating part 8 calculates the changing amount of 
the physical property values in the cases of taking the stress 
distribution into consideration and not taking the stress 
distribution into consideration (step S5). For example, 
with regard to the mobility of the carrier, the physical 
property value changing amount calculating part 8 calculates 
the changing amount ju( a) / jU of the mobility U( &) in the case 
of taking the stress distribution into consideration relative 
to the mobility jJL in the case of not taking the stress 
distribution into consideration. With regard to the 
diffusion coefficient, it calculates the changing amount 
D(CT)/D of the diffusion coefficient D(c7) in the case of 
taking the stress distribution relative to the diffusion 
coefficient D in the case of not taking the stress 
distribution into consideration. With regard to the viscous 
coefficient, it calculates the changing amount V{0)/V of 
the viscous coefficient ?? ( O) in the case of taking the stress 
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distribution into consideration relative to the viscous 
coefficient 7? in the case of not taking the stress 
distribution into consideration. 

Subsequently, the visualizing display control part 9 
5 displays the calculated result of the step S5 on the display 
apparatus 4 ( step S4 ) . Fig . 3B is a diagram showing an example 
of the changing amount ju ( O ) / ju of the mobility of the 
electrons. As shown in Fig. 3B, the mobility of the electron 
becomes low by 1-2 % in the surface region due to an influence 

10 of the compression stress occurring in the silicon substrate 
right under the gate. 

Thus, the present embodiment visually displays the 
changing amount of the physical property value in the cases 
of taking the stress distribution into consideration and not 

15 taking the stress distribution into consideration. 
Accordingly, it is possible to quantitatively and visually 
estimate how a diffusion length changes when the physical 
property value is the diffusion coefficient, how a relaxation 
time changes when the physical property value is the viscous 

20 coefficient, or how a current value when the physical property 
value is the carrier mobility, due to the stress. 
Furthermore, it is possible to easily grasp which part in 
the semiconductor device is affected by each of the physical 
property values . 

25 In the above-mentioned embodiment, although the 

example of calculating the physical property value obtained 
by taking the stress into consideration has been described, 
the physical quantity to be taken into consideration is not 
limited to the stress. For example, the physical property 

30 value such as the diffusion coefficient may be calculated 
by taking concentrations of the electrons and the holes into 
consideration. 

(Second Embodiment) 

A second embodiment has a feature that the changing 
35 amount of the physical property value is displayed on the 
distribution of the physical distribution in a superimposing 
manner. 
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Because the second embodiment has a block configuration 
similar to that of Fig. 1, a detailed description will be 
omitted hereinafter. 

Fig. 4 is a flow chart showing processes of the second 
embodiment of a simulation apparatus according to the present 
invention. The steps S11-S15 of Fig. 4 are the same as the 
steps S1-S5 of Fig. 2. After the process of the step S15 has 
been finished, the visualizing display control part 9 
displays the changing amount such as PL(G)I p. or D(<J)/D of 
the physical property value on the distribution of the 
physical amount (step S16) in the superimposing manner. 

Fig. 5B is a diagram showing the process result of the 
step S16 of Fig. 4. Fig. 5B shows an example of selecting 
the diffusion coefficient as the physical property value, 
and a boron concentration as the physical amount. Fig. 5A 
shows the stress distribution in the silicon substrate right 
under the gate, similarly to Fig. 3A. The stress of Fig. 5A 
expresses a hydrostatic pressure obtained as an average value 
between the stress in a crosswise direction and the stress 
in a lengthwise direction. 

As shown in Fig. 5B, in the step S16 of Fig. 4, the 
changing amount of the diffusion coefficients of the boron 
(the region right under of the gate of Fig. 5) in the cases 
of taking the stress distribution into consideration and not 
taking the stress distribution into consideration is 
displayed on the concentration distribution of the boron (a 
contour line in a right-and-left direction of the gate in 
Fig. 5) in the superimposing manner. 

As shown in Fig. 5, although the diffusion coefficient 
of the boron increases by 1-2 % due to the compression stress 
right under the gate, the boron does not exist right under 
the gate almost at all. Accordingly, it may be unnecessary 
to take the stress into consideration. 

Thus, because the present embodiment displays the 
changing amount of the physical property values in the cases 
of taking the stress into consideration and not taking the 
stress into consideration on the distribution of the physical 



quantity in the superimposing manner, it is possible to easily 
and accurately estimate whether or not the it is necessary 
to take the stress into consideration when the semiconductor 
device is developed. 

Furthermore, in the step SI 6 of Fig. 4, the changing 
amount of the viscous coefficients in the cases of taking 
the stress into consideration and not taking the stress into 
consideration may be displayed on the stress distribution 
in the superimposing manner. 

Furthermore, the mobility may be selected as the 
physical property value, the current path may be selected 
as the physical amount, and the changing amounts of the 
mobility in the cases of taking the stress distribution into 
consideration and not taking the stress distribution into 
consideration may be displayed on the distribution of the 
current path in the superimposing manner. 

(Third Embodiment) 

The third embodiment has a feature that a ratio of the 
changing amount of the physical property value at one 
fabrication condition relative to the changing amount of the 
physical property value at the other fabrication condition 
is calculated. 

Because the third embodiment has also a block 
configuration similar to that of Fig. 1, the detailed 
description will be omitted. 

Fig. 6 is a flow chart showing processes of the third 
embodiment of a simulation apparatus according to the third 
embodiment. First of all, the device simulation executing 
part 1 executes a simulation based on a first fabrication 
condition (step S21), and stores the simulated result (step 
S22) . In step S21, the changing amount between the physical 
property value in the case of taking the stress distribution 
into consideration and the physical property value in the 
case of not taking the stress distribution into consideration 
is calculated. 

Furthermore, the device simulation executing part 1 
executes a simulation based on a second fabrication condition 
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(step S23) , and stores the simulated result (step S24) . The 
changing amount of the physical property values in the cases 
of taking the stress distribution into consideration and not 
taking the stress distribution into consideration, based on 
5 the second fabrication condition. 

Which of the steps S21 and S22, and the steps S23 and 
S24 may be executed on ahead. 

Subsequently, the physical property value is selected 
(step S25), and the device shape obtained by the simulation 
10 at the first fabrication condition is displayed on the device 
shape obtained by the simulation at the second fabrication 
condition in the superimposing manner(step S26). 
JjT Subsequently, a ratio of the changing amount of the 

O physical property value calculated by the step S23 relative 

^L; 15 to the changing amount of the physical property value 

ifi calculated by the step S21 is calculated ( step S27). More 

't! specifically, ( jul2( G2) / jj.2) / { jul1{ <J2) / /ml) is calculated 

with regard to the carrier mobility, (D2( cr 2 ) /D2 )/( Dl ( G 
1 ) /Dl ) is calculated with regard to the diffusion coefficient, 
fy 20 and (772(cr2 )/772)/(7?l(cri)/77l) is calculated with regard 

y[l to the viscous coefficient. 

Si Subsequently, the calculated result of the step S27 is 

visually displayed on the display apparatus 4 (step S27). 
Fig. 7 is a diagram showing the process result of the 
25 step S27. Fig. 7A shows a stress distribution occurring in 
the silicon substrate of the bulk, and Fig. 7B is the stress 
distribution occurring in a SOI substrate. The stresses of 
Fig. 7A and Fig. 7B express a hydrostatic pressure obtained 
as an average value between the stress in the crosswise 
30 direction and the stress in the lengthwise direction. 

Fig. 7C shows the changing amount of the mobility of 
the electrons calculated from the stress distribution of Fig. 
7A and Fig. 7B. More specifically, Fig. 7C shows a ratio of 
the changing amount of the mobility of the electrons in the 
35 SOI substrate relative to the changing amount of the mobility 
of the electrons in the silicon substrate of the bulk. 

Because the stress value occurring in the SOI substrate 
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is higher than the stress value occurring in the silicon 
substrate, the mobility of the electron becomes low by 2-4 % 
in the SOI substrate relative to the silicon of the bulk. 

Thus, in the third embodiment, it is possible to compare 
the changing amounts of the physical property due to the 
stress with each other, with regard to a plurality of the 
semiconductor device experimentally fabricated on the SOI 
substrate, for example, it is possible to calculate the 
changing amount due to the stress of the semiconductor device 
experimentally fabricated on the SOI substrate. Therefore, 
it is possible to exactly grasp how the stress is affected 
when the fabrication processes are changed. Accordingly, it 
is possible to select the proper fabrication process, thereby 
reducing the fabrication cost and shortening the development 
period . 

In each of the above-mentioned embodiments, specified 
forms for visualization are not limited to the illustrated 
one. 

Furthermore, the simulation apparatus described in 
each of the above-mentioned embodiments may be composed of 
a hardware or a software. When composed of the software, a 
program for realizing functions of the simulation apparatus 
is contained in recording mediums such as a floppy disk 11 
or a CD-ROM 12 shown in Fig. 8, and may be read by the computer 
13 to carry out it. The recording medium is not limited to 
portable mediums such as a magnetic disk or an optical disk. 
Fixed type recording mediums such as a hard drive apparatus 
or a semiconductor memory are also applicable to the present 
invention. 

Furthermore, the program for realizing functions of the 
simulation apparatus may be delivered via a communication 
line (including a wireless communication line) such as an 
Internet. Furthermore, the program may be transmitted via 
a wired line or a wireless line such as the Internet at the 
encrypted, modulated, and compressed state, or by being 
contained in the recording medium. 
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